SpatialDDLS: an R package to deconvolute spatial transcriptomics data using neural networks.

SUMMARY: Spatial transcriptomics has changed our way to study tissue structure and cellular organization. However, there are still limitations in its resolution, and most available platforms do not reach a single cell resolution. To address this issue, we introduce SpatialDDLS, a fast neural network-based algorithm for cell type deconvolution of spatial transcriptomics data. SpatialDDLS leverages single-cell RNA sequencing data to simulate mixed transcriptional profiles with predefined cellular composition, which are subsequently used to train a fully connected neural network to uncover cell type diversity within each spot. By comparing it with two state-of-the-art spatial deconvolution methods, we demonstrate that SpatialDDLS is an accurate and fast alternative to the available state-of-the art tools. AVAILABILITY AND IMPLEMENTATION: The R package SpatialDDLS is available via CRAN-The Comprehensive R Archive Network: https://CRAN.R-project.org/package=SpatialDDLS. A detailed manual of the main functionalities implemented in the package can be found at https://diegommcc.github.io/SpatialDDLS.

Caveolin-1 dolines form a distinct and rapid caveolae-independent mechanoadaptation system.

In response to different types and intensities of mechanical force, cells modulate their physical properties and adapt their plasma membrane (PM). Caveolae are PM nano-invaginations that contribute to mechanoadaptation, buffering tension changes. However, whether core caveolar proteins contribute to PM tension accommodation independently from the caveolar assembly is unknown. Here we provide experimental and computational evidence supporting that caveolin-1 confers deformability and mechanoprotection independently from caveolae, through modulation of PM curvature. Freeze-fracture electron microscopy reveals that caveolin-1 stabilizes non-caveolar invaginations-dolines-capable of responding to low-medium mechanical forces, impacting downstream mechanotransduction and conferring mechanoprotection to cells devoid of caveolae. Upon cavin-1/PTRF binding, doline size is restricted and membrane buffering is limited to relatively high forces, capable of flattening caveolae. Thus, caveolae and dolines constitute two distinct albeit complementary components of a buffering system that allows cells to adapt efficiently to a broad range of mechanical stimuli.

A deep learning-based tool for the automated detection and analysis of caveolae in transmission electron microscopy images.

Caveolae are nanoscopic and mechanosensitive invaginations of the plasma membrane, essential for adipocyte biology. Transmission electron microscopy (TEM) offers the highest resolution for caveolae visualization, but provides complicated images that are difficult to classify or segment using traditional automated algorithms such as threshold-based methods. As a result, the time-consuming tasks of localization and quantification of caveolae are currently performed manually. We used the Keras library in R to train a convolutional neural network with a total of 36,000 TEM image crops obtained from adipocytes previously annotated manually by an expert. The resulting model can differentiate caveolae from non-caveolae regions with a 97.44% accuracy. The predictions of this model are further processed to obtain caveolae central coordinate detection and cytoplasm boundary delimitation. The model correctly finds negligible caveolae predictions in images from caveolae depleted Cav1-/- adipocytes. In large reconstructions of adipocyte sections, model and human performances are comparable. We thus provide a new tool for accurate caveolae automated analysis that could speed up and assist in the characterization of the cellular mechanical response.

Chromatin dynamics through mouse preimplantation development revealed by single molecule localisation microscopy.

Most studies addressing chromatin behaviour during preimplantation development are based on biochemical assays that lack spatial and cell-specific information, crucial during early development. Here, we describe the changes in chromatin taking place at the transition from totipotency to lineage specification, by using direct stochastical optical reconstruction microscopy (dSTORM) in whole-mount embryos during the first stages of mouse development. Through the study of two post-translational modifications of Histone 3 related to active and repressed chromatin, H3K4me3 and H3K9me3 respectively, we obtained a time-course of chromatin states, showing spatial differences between cell types, related to their differentiation state. This analysis adds a new layer of information to previous biochemical studies and provides novel insight to current models of chromatin organisation during the first stages of development.

RXRs control serous macrophage neonatal expansion and identity and contribute to ovarian cancer progression.

Tissue-resident macrophages (TRMs) populate all tissues and play key roles in homeostasis, immunity and repair. TRMs express a molecular program that is mostly shaped by tissue cues. However, TRM identity and the mechanisms that maintain TRMs in tissues remain poorly understood. We recently found that serous-cavity TRMs (LPMs) are highly enriched in RXR transcripts and RXR-response elements. Here, we show that RXRs control mouse serous-macrophage identity by regulating chromatin accessibility and the transcriptional regulation of canonical macrophage genes. RXR deficiency impairs neonatal expansion of the LPM pool and reduces the survival of adult LPMs through excess lipid accumulation. We also find that peritoneal LPMs infiltrate early ovarian tumours and that RXR deletion diminishes LPM accumulation in tumours and strongly reduces ovarian tumour progression in mice. Our study reveals that RXR signalling controls the maintenance of the serous macrophage pool and that targeting peritoneal LPMs may improve ovarian cancer outcomes.

Instantaneous and non-destructive relative water content estimation from deep learning applied to resonant ultrasonic spectra of plant leaves.

BACKGROUND: Non-contact resonant ultrasound spectroscopy (NC-RUS) has been proven as a reliable technique for the dynamic determination of leaf water status. It has been already tested in more than 50 plant species. In parallel, relative water content (RWC) is highly used in the ecophysiological field to describe the degree of water saturation in plant leaves. Obtaining RWC implies a cumbersome and destructive process that can introduce artefacts and cannot be determined instantaneously. RESULTS: Here, we present a method for the estimation of RWC in plant leaves from non-contact resonant ultrasound spectroscopy (NC-RUS) data. This technique enables to collect transmission coefficient in a [0.15-1.6] MHz frequency range from plant leaves in a non-invasive, non-destructive and rapid way. Two different approaches for the proposed method are evaluated: convolutional neural networks (CNN) and random forest (RF). While CNN takes the entire ultrasonic spectra acquired from the leaves, RF only uses four relevant parameters resulted from the transmission coefficient data. Both methods were tested successfully in Viburnum tinus leaf samples with Pearson's correlations between 0.92 and 0.84. CONCLUSIONS: This study showed that the combination of NC-RUS technique with deep learning algorithms is a robust tool for the instantaneous, accurate and non-destructive determination of RWC in plant leaves.

A graph-cut approach for pulmonary artery-vein segmentation in noncontrast CT images.

Lung vessel segmentation has been widely explored by the biomedical image processing community; however, the differentiation of arterial from venous irrigation is still a challenge. Pulmonary artery-vein (AV) segmentation using computed tomography (CT) is growing in importance owing to its undeniable utility in multiple cardiopulmonary pathological states, especially those implying vascular remodelling, allowing the study of both flow systems separately. We present a new framework to approach the separation of tree-like structures using local information and a specifically designed graph-cut methodology that ensures connectivity as well as the spatial and directional consistency of the derived subtrees. This framework has been applied to the pulmonary AV classification using a random forest (RF) pre-classifier to exploit the local anatomical differences of arteries and veins. The evaluation of the system was performed using 192 bronchopulmonary segment phantoms, 48 anthropomorphic pulmonary CT phantoms, and 26 lungs from noncontrast CT images with precise voxel-based reference standards obtained by manually labelling the vessel trees. The experiments reveal a relevant improvement in the accuracy ( ∼ 20%) of the vessel particle classification with the proposed framework with respect to using only the pre-classification based on local information applied to the whole area of the lung under study. The results demonstrated the accurate differentiation between arteries and veins in both clinical and synthetic cases, specifically when the image quality can guarantee a good airway segmentation, which opens a huge range of possibilities in the clinical study of cardiopulmonary diseases.

Tox_(R)CNN: Deep learning-based nuclei profiling tool for drug toxicity screening.

Toxicity is an important factor in failed drug development, and its efficient identification and prediction is a major challenge in drug discovery. We have explored the potential of microscopy images of fluorescently labeled nuclei for the prediction of toxicity based on nucleus pattern recognition. Deep learning algorithms obtain abstract representations of images through an automated process, allowing them to efficiently classify complex patterns, and have become the state-of-the art in machine learning for computer vision. Here, deep convolutional neural networks (CNN) were trained to predict toxicity from images of DAPI-stained cells pre-treated with a set of drugs with differing toxicity mechanisms. Different cropping strategies were used for training CNN models, the nuclei-cropping-based Tox_CNN model outperformed other models classifying cells according to health status. Tox_CNN allowed automated extraction of feature maps that clustered compounds according to mechanism of action. Moreover, fully automated region-based CNNs (RCNN) were implemented to detect and classify nuclei, providing per-cell toxicity prediction from raw screening images. We validated both Tox_(R)CNN models for detection of pre-lethal toxicity from nuclei images, which proved to be more sensitive and have broader specificity than established toxicity readouts. These models predicted toxicity of drugs with mechanisms of action other than those they had been trained for and were successfully transferred to other cell assays. The Tox_(R)CNN models thus provide robust, sensitive, and cost-effective tools for in vitro screening of drug-induced toxicity. These models can be adopted for compound prioritization in drug screening campaigns, and could thereby increase the efficiency of drug discovery.

Von Hippel-Lindau Protein Is Required for Optimal Alveolar Macrophage Terminal Differentiation, Self-Renewal, and Function.

The rapid transit from hypoxia to normoxia in the lung that follows the first breath in newborn mice coincides with alveolar macrophage (AM) differentiation. However, whether sensing of oxygen affects AM maturation and function has not been previously explored. We have generated mice whose AMs show a deficient ability to sense oxygen after birth by deleting Vhl, a negative regulator of HIF transcription factors, in the CD11c compartment (CD11cΔVhl mice). VHL-deficient AMs show an immature-like phenotype and an impaired self-renewal capacity in vivo that persists upon culture ex vivo. VHL-deficient phenotype is intrinsic in AMs derived from monocyte precursors in mixed bone marrow chimeras. Moreover, unlike control Vhlfl/fl, AMs from CD11cΔVhl mice do not reverse pulmonary alveolar proteinosis when transplanted into Csf2rb-/- mice, demonstrating that VHL contributes to AM-mediated surfactant clearance. Thus, our results suggest that optimal AM terminal differentiation, self-renewal, and homeostatic function requires their intact oxygen-sensing capacity.

Pulmonary Artery-Vein Classification in CT Images Using Deep Learning.

Recent studies show that pulmonary vascular diseases may specifically affect arteries or veins through different physiologic mechanisms. To detect changes in the two vascular trees, physicians manually analyze the chest computed tomography (CT) image of the patients in search of abnormalities. This process is time consuming, difficult to standardize, and thus not feasible for large clinical studies or useful in real-world clinical decision making. Therefore, automatic separation of arteries and veins in CT images is becoming of great interest, as it may help physicians to accurately diagnose pathological conditions. In this paper, we present a novel, fully automatic approach to classify vessels from chest CT images into arteries and veins. The algorithm follows three main steps: first, a scale-space particles segmentation to isolate vessels; then a 3-D convolutional neural network (CNN) to obtain a first classification of vessels; finally, graph-cuts' optimization to refine the results. To justify the usage of the proposed CNN architecture, we compared different 2-D and 3-D CNNs that may use local information from bronchus- and vessel-enhanced images provided to the network with different strategies. We also compared the proposed CNN approach with a random forests (RFs) classifier. The methodology was trained and evaluated on the superior and inferior lobes of the right lung of 18 clinical cases with noncontrast chest CT scans, in comparison with manual classification. The proposed algorithm achieves an overall accuracy of 94%, which is higher than the accuracy obtained using other CNN architectures and RF. Our method was also validated with contrast-enhanced CT scans of patients with chronic thromboembolic pulmonary hypertension to demonstrate that our model generalizes well to contrast-enhanced modalities. The proposed method outperforms state-of-the-art methods, paving the way for future use of 3-D CNN for artery/vein classification in CT images.

Flow Cytometry Data Preparation Guidelines for Improved Automated Phenotypic Analysis.

Advances in flow cytometry (FCM) increasingly demand adoption of computational analysis tools to tackle the ever-growing data dimensionality. In this study, we tested different data input modes to evaluate how cytometry acquisition configuration and data compensation procedures affect the performance of unsupervised phenotyping tools. An analysis workflow was set up and tested for the detection of changes in reference bead subsets and in a rare subpopulation of murine lymph node CD103+ dendritic cells acquired by conventional or spectral cytometry. Raw spectral data or pseudospectral data acquired with the full set of available detectors by conventional cytometry consistently outperformed datasets acquired and compensated according to FCM standards. Our results thus challenge the paradigm of one-fluorochrome/one-parameter acquisition in FCM for unsupervised cluster-based analysis. Instead, we propose to configure instrument acquisition to use all available fluorescence detectors and to avoid integration and compensation procedures, thereby using raw spectral or pseudospectral data for improved automated phenotypic analysis.

Pluripotency Surveillance by Myc-Driven Competitive Elimination of Differentiating Cells.

The mammalian epiblast is formed by pluripotent cells able to differentiate into all tissues of the new individual. In their progression to differentiation, epiblast cells and their in vitro counterparts, embryonic stem cells (ESCs), transit from naive pluripotency through a differentiation-primed pluripotent state. During these events, epiblast cells and ESCs are prone to death, driven by competition between Myc-high cells (winners) and Myc-low cells (losers). Using live tracking of Myc levels, we show that Myc-high ESCs approach the naive pluripotency state, whereas Myc-low ESCs are closer to the differentiation-primed state. In ESC colonies, naive cells eliminate differentiating cells by cell competition, which is determined by a limitation in the time losers are able to survive persistent contact with winners. In the mouse embryo, cell competition promotes pluripotency maintenance by elimination of primed lineages before gastrulation. The mechanism described here is relevant to mammalian embryo development and induced pluripotency.

ESC-Track: A computer workflow for 4-D segmentation, tracking, lineage tracing and dynamic context analysis of ESCs.

Embryonic stem cells (ESCs) can be established as permanent cell lines, and their potential to differentiate into adult tissues has led to widespread use for studying the mechanisms and dynamics of stem cell differentiation and exploring strategies for tissue repair. Imaging live ESCs during development is now feasible due to advances in optical imaging and engineering of genetically encoded fluorescent reporters; however, a major limitation is the low spatio-temporal resolution of long-term 3-D imaging required for generational and neighboring reconstructions. Here, we present the ESC-Track (ESC-T) workflow, which includes an automated cell and nuclear segmentation and tracking tool for 4-D (3-D + time) confocal image data sets as well as a manual editing tool for visual inspection and error correction. ESC-T automatically identifies cell divisions and membrane contacts for lineage tree and neighborhood reconstruction and computes quantitative features from individual cell entities, enabling analysis of fluorescence signal dynamics and tracking of cell morphology and motion. We use ESC-T to examine Myc intensity fluctuations in the context of mouse ESC (mESC) lineage and neighborhood relationships. ESC-T is a powerful tool for evaluation of the genealogical and microenvironmental cues that maintain ESC fitness.

Implementation and Performance of Automated Software for Computing Right-to-Left Ventricular Diameter Ratio From Computed Tomography Pulmonary Angiography Images.

OBJECTIVE: The aim of this study was to prospectively test the performance and potential for clinical integration of software that automatically calculates the right-to-left ventricular (RV/LV) diameter ratio from computed tomography pulmonary angiography images. METHODS: Using 115 computed tomography pulmonary angiography images that were positive for acute pulmonary embolism, we prospectively evaluated RV/LV ratio measurements that were obtained as follows: (1) completely manual measurement (reference standard), (2) completely automated measurement using the software, and (3 and 4) using a customized software interface that allowed 2 independent radiologists to manually adjust the automatically positioned calipers. RESULTS: Automated measurements underestimated (P < 0.001) the reference standard (1.09 [0.25] vs1.03 [0.35]). With manual correction of the automatically positioned calipers, the mean ratio became closer to the reference standard (1.06 [0.29] by read 1 and 1.07 [0.30] by read 2), and the correlation improved (r = 0.675 to 0.872 and 0.887). The mean time required for manual adjustment (37 [20] seconds) was significantly less than the time required to perform measurements entirely manually (100 [23] seconds). CONCLUSIONS: Automated CT RV/LV diameter ratio software shows promise for integration into the clinical workflow for patients with acute pulmonary embolism.

Automatic Synthesis of Anthropomorphic Pulmonary CT Phantoms.

The great density and structural complexity of pulmonary vessels and airways impose limitations on the generation of accurate reference standards, which are critical in training and in the validation of image processing methods for features such as pulmonary vessel segmentation or artery-vein (AV) separations. The design of synthetic computed tomography (CT) images of the lung could overcome these difficulties by providing a database of pseudorealistic cases in a constrained and controlled scenario where each part of the image is differentiated unequivocally. This work demonstrates a complete framework to generate computational anthropomorphic CT phantoms of the human lung automatically. Starting from biological and image-based knowledge about the topology and relationships between structures, the system is able to generate synthetic pulmonary arteries, veins, and airways using iterative growth methods that can be merged into a final simulated lung with realistic features. A dataset of 24 labeled anthropomorphic pulmonary CT phantoms were synthesized with the proposed system. Visual examination and quantitative measurements of intensity distributions, dispersion of structures and relationships between pulmonary air and blood flow systems show good correspondence between real and synthetic lungs (p > 0.05 with low Cohen's d effect size and AUC values), supporting the potentiality of the tool and the usefulness of the generated phantoms in the biomedical image processing field.

Automated axial right ventricle to left ventricle diameter ratio computation in computed tomography pulmonary angiography.

BACKGROUND AND PURPOSE: Right Ventricular to Left Ventricular (RV/LV) diameter ratio has been shown to be a prognostic biomarker for patients suffering from acute Pulmonary Embolism (PE). While Computed Tomography Pulmonary Angiography (CTPA) images used to confirm a clinical suspicion of PE do include information of the heart, a numerical RV/LV diameter ratio is not universally reported, likely because of lack in training, inter-reader variability in the measurements, and additional effort by the radiologist. This study designs and validates a completely automated Computer Aided Detection (CAD) system to compute the axial RV/LV diameter ratio from CTPA images so that the RV/LV diameter ratio can be a more objective metric that is consistently reported in patients for whom CTPA diagnoses PE. MATERIALS AND METHODS: The CAD system was designed specifically for RV/LV measurements. The system was tested in 198 consecutive CTPA patients with acute PE. Its accuracy was evaluated using reference standard RV/LV radiologist measurements and its prognostic value was established for 30-day PE-specific mortality and a composite outcome of 30-day PE-specific mortality or the need for intensive therapies. The study was Institutional Review Board (IRB) approved and HIPAA compliant. RESULTS: The CAD system analyzed correctly 92.4% (183/198) of CTPA studies. The mean difference between automated and manually computed axial RV/LV ratios was 0.03±0.22. The correlation between the RV/LV diameter ratio obtained by the CAD system and that obtained by the radiologist was high (r=0.81). Compared to the radiologist, the CAD system equally achieved high accuracy for the composite outcome, with areas under the receiver operating characteristic curves of 0.75 vs. 0.78. Similar results were found for 30-days PE-specific mortality, with areas under the curve of 0.72 vs. 0.75. CONCLUSIONS: An automated CAD system for determining the CT derived RV/LV diameter ratio in patients with acute PE has high accuracy when compared to manual measurements and similar prognostic significance for two clinical outcomes.

Right ventricle segmentation from cardiac MRI: a collation study.

Magnetic Resonance Imaging (MRI), a reference examination for cardiac morphology and function in humans, allows to image the cardiac right ventricle (RV) with high spatial resolution. The segmentation of the RV is a difficult task due to the variable shape of the RV and its ill-defined borders in these images. The aim of this paper is to evaluate several RV segmentation algorithms on common data. More precisely, we report here the results of the Right Ventricle Segmentation Challenge (RVSC), concretized during the MICCAI'12 Conference with an on-site competition. Seven automated and semi-automated methods have been considered, along them three atlas-based methods, two prior based methods, and two prior-free, image-driven methods that make use of cardiac motion. The obtained contours were compared against a manual tracing by an expert cardiac radiologist, taken as a reference, using Dice metric and Hausdorff distance. We herein describe the cardiac data composed of 48 patients, the evaluation protocol and the results. Best results show that an average 80% Dice accuracy and a 1cm Hausdorff distance can be expected from semi-automated algorithms for this challenging task on the datasets, and that an automated algorithm can reach similar performance, at the expense of a high computational burden. Data are now publicly available and the website remains open for new submissions (http://www.litislab.eu/rvsc/).

Comparing algorithms for automated vessel segmentation in computed tomography scans of the lung: the VESSEL12 study.

The VESSEL12 (VESsel SEgmentation in the Lung) challenge objectively compares the performance of different algorithms to identify vessels in thoracic computed tomography (CT) scans. Vessel segmentation is fundamental in computer aided processing of data generated by 3D imaging modalities. As manual vessel segmentation is prohibitively time consuming, any real world application requires some form of automation. Several approaches exist for automated vessel segmentation, but judging their relative merits is difficult due to a lack of standardized evaluation. We present an annotated reference dataset containing 20 CT scans and propose nine categories to perform a comprehensive evaluation of vessel segmentation algorithms from both academia and industry. Twenty algorithms participated in the VESSEL12 challenge, held at International Symposium on Biomedical Imaging (ISBI) 2012. All results have been published at the VESSEL12 website http://vessel12.grand-challenge.org. The challenge remains ongoing and open to new participants. Our three contributions are: (1) an annotated reference dataset available online for evaluation of new algorithms; (2) a quantitative scoring system for objective comparison of algorithms; and (3) performance analysis of the strengths and weaknesses of the various vessel segmentation methods in the presence of various lung diseases.

A benchmark for comparison of cell tracking algorithms.

MOTIVATION: Automatic tracking of cells in multidimensional time-lapse fluorescence microscopy is an important task in many biomedical applications. A novel framework for objective evaluation of cell tracking algorithms has been established under the auspices of the IEEE International Symposium on Biomedical Imaging 2013 Cell Tracking Challenge. In this article, we present the logistics, datasets, methods and results of the challenge and lay down the principles for future uses of this benchmark. RESULTS: The main contributions of the challenge include the creation of a comprehensive video dataset repository and the definition of objective measures for comparison and ranking of the algorithms. With this benchmark, six algorithms covering a variety of segmentation and tracking paradigms have been compared and ranked based on their performance on both synthetic and real datasets. Given the diversity of the datasets, we do not declare a single winner of the challenge. Instead, we present and discuss the results for each individual dataset separately. AVAILABILITY AND IMPLEMENTATION: The challenge Web site (http://www.codesolorzano.com/celltrackingchallenge) provides access to the training and competition datasets, along with the ground truth of the training videos. It also provides access to Windows and Linux executable files of the evaluation software and most of the algorithms that competed in the challenge.

Organ-focused mutual information for nonrigid multimodal registration of liver CT and Gd-EOB-DTPA-enhanced MRI.

Accurate detection of liver lesions is of great importance in hepatic surgery planning. Recent studies have shown that the detection rate of liver lesions is significantly higher in gadoxetic acid-enhanced magnetic resonance imaging (Gd-EOB-DTPA-enhanced MRI) than in contrast-enhanced portal-phase computed tomography (CT); however, the latter remains essential because of its high specificity, good performance in estimating liver volumes and better vessel visibility. To characterize liver lesions using both the above image modalities, we propose a multimodal nonrigid registration framework using organ-focused mutual information (OF-MI). This proposal tries to improve mutual information (MI) based registration by adding spatial information, benefiting from the availability of expert liver segmentation in clinical protocols. The incorporation of an additional information channel containing liver segmentation information was studied. A dataset of real clinical images and simulated images was used in the validation process. A Gd-EOB-DTPA-enhanced MRI simulation framework is presented. To evaluate results, warping index errors were calculated for the simulated data, and landmark-based and surface-based errors were calculated for the real data. An improvement of the registration accuracy for OF-MI as compared with MI was found for both simulated and real datasets. Statistical significance of the difference was tested and confirmed in the simulated dataset (p<0.01).